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With the reduction of funding for sea ice reconnaissance flights, the 
National/Naval Ice Center needs to capitalize on the improvements in satellite 
technology. Imaging sensors such as AVHRR, DMSP/OLS, SSM/I and RADARSAT are 
used to detect the presence of sea ice, but with the exception of SSM/I and RADARSAT, 
clouds are a major obstacle to viewing the surface. With NASA’s development of the 
Moderate-resolution Imaging Spectroradiometer (MODIS) and MODIS Airborne 
Simulator (MAS), there is finally a sensor capable of using multi-spectral techniques to 
detect the presence of clouds. 

A group at the Space Science and Engineering Center (SSEC), University of 
Wisconsin - Madison lead by Dr. Steve Ackerman has developed a cloud mask for 
MAS/MODIS. The technique determines a level of confidence that a given pixel is clear 
based on a series of multi-spectral tests. By combining the confidence level from all 
tests, it is possible to detect the presence of clouds at different altitudes in the 
atmosphere. Threshold optimizations are described in this thesis for the T B (ll|J.m) and 
T B (3.9pm) - T B (ll!im) tests from Ackerman et al. (1997). In addition, the T B (ll|im) - 
T B (12pm) test is removed. These modifications are based on daytime analysis of several 
MAS cases and a limited number of MODIS cases. 

Subjective analysis shows the modifications greatly improve the detection of 
clouds over cold polar oceans where sub-pixel ice may be present or water temperatures 
might falsely indicate clouds. The number of Cloudy pixels (<0.66 clear confidence 
level) for a given scene was increased 12.1% on average for MAS cases. The NPS cloud 
mask also classified two times more Probably Clear and Undecided pixels than the 
original mask due to greater sensitivity to thin, small clouds. 
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INTRODUCTION 



A. MOTIVATION 

With the reduction of funding for sea ice reconnaissance flights, manned polar 
observation stations, and polar in situ measurements, the National/Naval Ice Center 
(NIC) needs to capitalize on the latest advances in remote sensing. Clouds can be 
characterized by either their brightness temperature or reflectance characteristics. 
Previous cloud masks used simple combinations of visible and IR wavelengths to detect 
the presence of clouds. Cloud detection is difficult in polar ocean regions where sea ice 
may have the same brightness temperature and/or reflectance of the surrounding clouds. 
Detection becomes further complicated when clouds are thin along their edges and are 
almost indistinguishable from the sea ice. 

The Moderate Resolution Imaging Spectroradiometer (MODIS) and the MODIS 
Airborne Simulator (MAS) provide an opportunity for multi-spectral cloud detection 
approaches using many more wavelengths than previously available. Multi-spectral 
cloud detection will take advantage of the many differences in the spectral properties of 
sea ice and clouds. As a result, future automated sea ice analysis algorithms will more 
correctly identify sea ice properties with minimal cloud contamination. 

B. HISTORY 

Operational analysis of sea ice using remote sensing assets date back to the 
1970’s. Analysts manually analyzed sea ice in imagery with grease pencils and 
eventually transferred the analyses to paper charts. In the early 1990’s, on-screen sea ice 
analysis became possible due to increased computer power and products were 
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disseminated via the World Wide Web. The future objective of operational sea ice 
analysis is an automated process with minimal manual intervention. Automated analysis 
requires a cloud mask that is sensitive to the unique conditions present in polar ocean 
regions. 

Satellite sensors are available that can penetrate clouds to derive surface 
information without the need of a cloud mask. One is a passive microwave sensor called 
Special Sensor Microwave Imager (SSM/I) on the Defense Meteorological Satellite 
Program (DMSP). SSM/I techniques usually use 85GHz imagery to identify sea ice and 
cloud liquid water content. SSM/I is not the perfect ice detection instrument. Microwave 
techniques currently provide low-resolution observations (ca. 25km) and have difficulty 
identifying sea ice concentrations when water pools on top of the sea ice. A second 
active technique uses the Synthetic Aperture Radar (SAR). For example RADARSAT 
uses a 5.3GHz radar to penetrate clouds and view the surface. SAR is an ideal tool for 
operational sea ice analysis but is restrictive due to the cost per scene to the U.S. 
Government. 

C. OBJECTIVE 

The objective of this thesis is to optimize a cloud mask for MODIS that is 
optimized for the polar ocean region and used by NIC and other polar scientists. The 
work described here is based on a cloud mask that was originally developed by the Space 
Science and Engineering Center at the University of Wisconsin - Madison. The 
optimized cloud mask will assist in the operational sea ice analysis by identifying clouds 
that might be mistaken as sea ice or their shadows mistaken as leads. The ultimate goal is 

the initialization of an automated sea ice analysis algorithm. 
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D. ORGANIZATION 



Chapter II outlines background information about the MAS/MODIS cloud mask. 
Chapter III describes the methodology for optimizing the MAS/MODIS Cloud mask. 
Chapter IV shows MAS examples of the NPS cloud mask adjustments. Finally, Chapter 
V shows a MODIS example of the NPS cloud mask and Chapter VI outlines the 
conclusions of this thesis and offers recommendations for further studies. 
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II. BACKGROUND INFORMATION 



A. PLATFORMS AND SENSORS 

Two platform/sensor packages are employed by the cloud mask algorithm. First 
is the ER-2 and MAS research sensor, which is used to evaluate channel properties. 
Second is the Terra Satellite and MODIS operational space-borne sensor. 

1. ER-2 and MAS 

MAS is a modified Daedalus Wildfire scanning spectrometer that is flown on a 
NASA ER-2 aircraft. It has been used in many experiments including WINCE (January 
- February 1997) and FIRE/ ACE (May - June 1998). It provides 50 band, 50m 
resolution from 20km flight altitude in 50 spectral bands range from 0.523pm to 
14.383pm wavelengths (Table 2-1). Band spectral ranges can change depending on the 
experiment. 

2. Terra Satellite and MODIS 

The Terra satellite, formally AM-1, is the flagship of the Earth Observing System 
(EOS). EOS is a major program under the National Aeronautics and Space 
Administration’s (NASA) Earth Science Enterprise (ESE) with a primary mission to 
observe the Earth’s surface. Terra was launched from the Western Test Range, 
Vandenburg Air Force Base, California on December 18, 1999 at 13:57 EST (Herring, 
2000). It is in a 10:30am local equator crossing, descending, Sun synchronous orbit in 
attempt to maximize the view of the surface while minimizing cloud cover. The orbit 
will be adjusted to ensure a complete surface coverage every 16 days with a life span of 6 
years (NASA, 1999a). 
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Table 2. 1 . MAS Instrument/Platform Specifications 



Platform 


NASA ER-2 
Aircraft 


Swath Width 


73.25 km or 22.9 miles 
(at 20 km altitude) 


Ground Speed 


400 kts (206 m/s) 


Total Field of 
View 


85.92° 


Pixel Spatial 
Resolution 


50 m (at 20 km) 


Instantaneous 
Field of view 


2.5 milliradians 


Pixel per Scan Line 


716 (roll corrected) 


Roll Correction 


+/- 3.5° 


Scan Rate 


6.25 scans/second 


Data Channels 


50 channels via digital 
recorded 


Spectral Bands 


Band 


Bandwidth 


Band 


Bandwidth 


1 


0.523 - 0.569 


26 


2.875 - 3.039 


2 


0.627 - 0.68 


27 


3.031-3.195 


3 


0.679-0.721 


28 


3.202-3.358 


4 


0.72-0.762 


29 


3.338-3.507 


5 


0.761 -0.802 


30 


3.511-3.668 


6 


0.802 - 0.844 


31 


3.669-3.819 


7 


0.844 - 0.886 


32 


3.819-3.986 


8 


0.888 - 0.922 


33 


3.979-4.132 


9 


0.921 -0.967 


34 


4.126-4.288 


10 


1.595 - 1.65 


35 


4.288 - 4.435 


11 


1.652-1.705 


36 


4.439-4.594 


12 


1.706-1.757 


37 


4.585-4.748 


13 


1.759- 1.812 


38 


4.745 - 4.901 


14 


1.813-1.863 


39 


4.897-5.047 


15 


1.863- 1.914 


40 


5.046 -5.203 


16 


1.914 -.1963 


41 


5.205 - 5.362 


17 


1.964-2.013 


42 


8.386-8.822 


18 


2.014-2.064 


43 


9.484-10.101 


19 


2.066-2.113 


44 


10.3 - 10.794 


20 


2.114-2.162 


45 


10.746-11.289 


21 


2.163-2.21 


46 


11.736-12.185 


22 


2.211-2.261 


47 


12.651 - 13.47 


23 


2.262-2.31 


48 


12.998-13.47 


24 


2.312-2.36 


49 


13.421 - 14.021 


25 


2.361-2.41 


50 


13.96- 14.383 



6 




One of the five instruments aboard TERRA is MODIS. It has 36 bands ranging 
from 0.4jim to 14.4{im wavelength. The bands are designed to observe land and ocean 
temperatures, land surface coverage, aerosols, water vapor, temperature profiles, and 
fires. Table 2-2 lists the various spectral bands, their associated spatial resolutions and 
primary use. 

B. MAS/MODIS ARCTIC CLOUD MASK 

The MAS/MODIS Cloud Mask is described in the algorithm theoretical basic 
document (ATBD) (Ackerman et al., 1997). The physical rationale and detail 
specifications of the cloud mask are documented in the ATBD. This thesis will only 
describe the important factors that effect cloud detection in polar ocean regions. 

Figure 2.1 shows MAS and MODIS bands overlaid on the transmittance 
properties for ozone, water vapor and carbon dioxide for a mid-latitude summer 
atmosphere calculated by MODTRAN 3.5 (Acharya et al., 1993). In general, cloud and 
surface observations are made in spectral regions with minimal atmospheric absorption. 
However, absorption bands can be used to limit the radiance from lower altitudes. This 
can help identify upper-level clouds. 

1. General overview of the Cloud Mask 

It is evident that the MAS and MODIS sensors have much in common in the 
exploitation of transmission and absorption features. In some cases they used different 
wavelengths to achieve the same result. As an example, MAS uses a 1.88fim band to 
exploit low and mid level strong water vapor absorption to highlight the presence of high, 
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Table 2.2. N 


10DIS Characteristics (NASA, 1999b) 


Primary Use 


Band 


Bandwidth 


Spatial Resolution 
(Nadir) 


Land/Cloud 

Boundaries 


1 


0.620- 0.670 urn 


250 m 


2 


0.841 -0.876 pm 


250 m 


Land/Cloud 

Properties 


3 


0.459 -0.479 pm 


500 m 


4 


0.545 -0.565 pm 


500 m 


5 


1.230- 1.250 pm 


500 m 


6 


1.628 - 1.652 pm 


500 m 


7 


2.105 -2.155 pm 


500 m 


Ocean Color/ 

Phytoplankton/ 

Biogeochemistry 


8 


0.405-0.420 pm 


1000 m 


9 


0.438 - 0.448 pm 


1000 m 


10 


0.483 - 0.493 pm 


1000 m 


11 


0.526 -0.536 pm 


1000 m 


12 


0.546 -0.556 pm 


1000 m 


13 


0.662 - 0.672 pm 


1000 m 


14 


0.673 - 0.683 pm 


1000 m 


15 


0.743 - 0.753 pm 


1000 m 


16 


0.862 -0.877 pm 


1000 m 


Atmospheric 
Water Vapor 


17 


0.890 -0.920 pm 


1000 m 


18 


0.931-0.941 pm 


1000 m 


19 


0.915-0.965 pm 


1000 m 


Surface/Cloud 

Temperature 


20 


3.660-3.840 pm 


1000 m 


21 


3.929 -3.989 pm 


1000 m 


22 


3.929 -3.989 pm 


1000 m 


23 


4.020 -4.080 pm 


1000 m 


Atmospheric 

Temperature 


24 


4.433 - 4.498 pm 


1000 m 


25 


4.482 - 4.549 pm 


1000 m 


Cirrus Clouds 


26 


1.360- 1.390 pm 


1000 m 


Water Vapor 


27 


6.535 - 6.895 pm 


1000 m 


28 


7.175 -7.475 pm 


1000 m 


29 


8.400 -8.700 pm 


1000 m 


Ozone 


30 


9.580- 9.880 pm 


1000 m 


Surface/Cloud 

Temperature 


31 


10.780- 11.280 pm 


1000 m 


32 


11.770 -12.270 pm 


1000 m 


Cloud Top 
Altitude 


33 


13.185 - 13.485 pm 


1000 m 


34 


13.485 - 13.785 pm 


1000 m 


35 


13.785 - 14.085 pm 


1000 m 


36 


. 14.085 -14.385 pm 


1000 m 
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Figure 2.1. MAS/MODIS bands - frequency relationship where these instruments exploit the various atmospheric properties to 
obtain information about the atmosphere. Overlaid are the transmittance properties for ozone, water vapor and 
carbon dioxide 



thin cirrus clouds. To accomplish the same task, MODIS uses the water vapor absorption 
band at 1.38p.m. 

There are five major groups of multi spectral tests in the mask. Each group of 
tests are designed to detect common cloud features ranging from low clouds to high 
cirrus clouds. The number of tests within a group and the number of groups used in a 
particular case vary depending on the time of day, the ecosystem, presence of snow or sea 
ice, amount of perceptible water, and several others. Once an individual pixel’s 
characteristics have been described (e.g. daytime, ocean, polar, ice), tests within various 
groups will be preformed along with specific thresholds for that pixel’s characteristics. 
All tests within a group and all groups have the same weight as any other test or group 
used for a pixel. 

a. Calculation of Confidence 

Ackerman et al. (1997) developed a simple linear scaling of thresholds. 
This technique allows for a degree of uncertainty to be included in each test. Figure 2.2 
depicts the linear function used by the cloud mask. There are three points; alpha (a) - 
any value below this point is considered cloudy, beta ((3) - this is the 50% clear 
confidence value (the pixel can be either clear or cloudy) and gamma (y) - any value 
above this point is considered clear. Each test within a group is given a value between 0 
and 1 depending on its threshold. A value of 0 is an indication of a high probability of 
cloud, while a value of 1 is an indication of a high probability of a clear condition. The 
lowest confidence value from the tests within a group becomes the value for the over all 
group. 
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The overall clear confidence value for the pixel, Q, is computed for the 



pixel as the N th root the product of group confidence values: 



Q=Jn s i- 

where Gi is the value for a group and N is the number of groups used in the test. This 
calculation is conservative in the estimate of a clear sky condition and will err on the 
cloudy side for any given pixel. For example, if a pixel were in a daytime, ocean, polar, 
sea ice location, the cloud mask would use four groups for the test. Confidence values 
for the four individual groups might be 0.82, 0.66, 1, and 1. Therefore the over all 
confidence would be [(0.82)(0.66)(1)(1)] /4 = 0.86 for that pixel. Table 2.3 shows Q 
values for the four confidence levels. The pixel in this example would be assigned as 
Undecided classification. 
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Table 2.3. Confidence Values 



Confidence 


Q 


Confident Clear 


Q > 0.99 


Probably Clear 


Q > 0.95 


l ndecidcd 


Q > 0.66 


Cloudy or Obstructed 


Q < 0.66 



2. Inputs 

The cloud mask requires some information about the scene before it can 
determine if a pixel is clear or not. These inputs determine the number of groups and the 
number of tests that will be preformed within each group. Some of the variables that are 
important for the polar ocean portion of the mask are day/night, ecosystem, and sea ice 
detection. 

a. Day/Night 

The day/night calculation is based on the solar zenith sun angle (0). When 
0 > 85° it is considered day time (Ackerman et ah, 1997). 

b. Ecosystem 

The ecosystem is determined by the pixel’s location within various index 

files and may be manually selected. They are: 

1 . Global USGS 1 km resolution land/sea tag 

2. Loveland 1 km resolution N. American ecosystem 

3. 10 minute resolution global Olson ecosystem 

4. Global Loveland 1 km resolution Olson ecosystem 

For the MAS cases, the Global Loveland 1 km resolution Olson ecosystem 
file was used to identify the ecosystem. The ecosystem is identified as either being 
desert, land, coast or ocean. 
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Sea Ice Detection 



c. 

There are challenges to overcome when it comes to the detection of sea 
ice. First, sea ice detection is difficult when there are low clouds in the scene due to their 
similar radiance and thermal signatures. Second, relatively thin sea ice detection is 
difficult because polar ocean waters might have similar radiance values. 

The MODIS cloud mask has two ways to detect the presence of sea ice. 
First, Special Sensor Microwave Imager (SSM/I) analysis at 25km resolution, available 
from either the National Center for Environmental Predictions (NCEP) or the National 
Snow and Ice Data Center (NSIDC), is used to identify the presence of sea ice. The 
inputs from NCEP and NSIDC are given equal weight in the mask for the determination 
of sea ice extent. Unfortunilty, SSM/I has difficulty in the summer when liquid water 
pools on top of sea ice and might incorrectly indicate an open water condition. 

Second, an internal means for idenfitying sea ice was developed by Hall et 
al., (1998) (herafter referred to the Hall Test) and implimented in the mask in an attempt 
to determine the pixel-scale extent in sea ice. In the Hall Test, when the brightness 
temperature at 11 Jim, Tb(1 1pm), is less than 280°K, the surface is assumed to be cold 
enough to have snow or sea ice present. The difference between snow or sea ice depends 
on whether the ecosystem is either land/coastal or ocean. Next, the Normalized 
Difference Snow Index (NDSI) test detects sea ice when its thickness is greater than 
10cm. This works because sea ice absorbs 1.62jim wavelength radiance. Therefore, 
when sea ice is present, visible reflectance, po.s5nm, will be greater than 1.62|im 
reflectance, pi. 62 nm- In addition, sea ice has a high reflectance at 0.88|im wavelength 
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(Po.88nm) and po.ssum can be used to identify its presence. Hall et al. ( 1 998) defines the 
NDSI as 

Po .55 jjm P\.62fiim 

When sea ice is present, Tb( 1 1pm) will be below 280°K and the NDSI will be computed. 
If the NDSI > 0.4 and po.88nm > 0- 1 then sea ice is possibly present. 

Finally, a test is needed to ensure that the assumed sea ice is not a high, 
cold cloud. Two bands are included in the Hall Test that are sensitive to high clouds. 
The atmosphere has a low tranmittance at 13pm due to CO 2 absorption and at 1.38pm 
due to water vapor absorption. In both cases, radiance from the lower troposphere is 
removed by absorption and only radiance emitted or reflected from the upper troposphere 
contributes to the radiance measured by the sensor. If T B (13pm) < 220°K and po^m > 
0.035, the Hall test assumes that high cloud is present in the pixel and is responsible for 
the false-positive assignment of sea ice. 

When applying the Hall Test, it is important to remember that polar ocean 
water has a very low albedo compared to that of the highly reflective arctic sea ice. 
Complications arise when dealing with “New Ice” (Frazil, Grease, Slush or Shuga) that 
has a thickness between 0 and 10 cm. There might not be enough contrast in reflectance 
(Po.ssum) to adequately separate the sea ice from water. In addition, as sea ice melts, 
water will pool on top of the floe and lower the albedo so the pixel will appear to be 
water when, in fact, it is sea ice. 
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For the operational MODIS cloud mask, sea ice is present only when 
either the NCEP or the NSIDC sea ice product and when the Hall Test detects the of sea 
ice. MAS only uses the Hall Test to detect the presence of sea ice. 

3. Groups and Test 

The five groups within the cloud mask are comprised of various tests. The 
number of groups and tests are determined by the input values. The groups that are 
important to the ocean polar regions are described below. 
a. Group 1 

The objective of Group 1 is to identify thick high clouds based on infrared 
Tb characteristics. This is accomplished in the ocean polar regions with Tb( 1 1 (lm) and 
T B (13.9jJ.m) wavelengths tests. The group has a potential difficulty with low clouds that 
have small thermal contrast between the cloud and the surface (Ackerman et al., 1997). 

(1) T b (1 l|im) Test. This is a simple IR cloud test that can obtain 
information about the surface or cloud top temperature. It is also important when dealing 
with the detection of snow covered sea ice as described above. In the open ocean when 
T e (l l|im) channel is expected to be less than 270°K ,the mask assumes the pixel is too 
cold to be open ocean and determines that it has failed the clear sky condition. In the 
polar ocean regions, where there might be subpixel sea ice present, there is a need to 
lower the T B (ll|im) threshold values. The magnitude of this modification will be 
discussed in the next chapter. In general, simple thresholds should vary with location, 
season, and time of day. 

(2) T B (13fim) Test. This test relies on the absorption 

characteristics of CO 2 . Clouds being sensed by a CO 2 band depend upon the atmospheric 
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attenuation in that band and the altitude of the cloud. The spectral band at 13.9pm 
provides good sensitivity to the relatively cold regions of the atmosphere. Only clouds 
above 500 hPa will have strong contributions to the radiance at the top of the atmosphere 
observed at 13.9pm. Lower regions of the atmosphere and the earth’s surface should 
provide negligible contributions to the radiance at 13.9 pm wavelength (Ackerman et al., 
1997). 

b. Group 2 

The objective of Group 2 is to detect thin clouds missed by Group 1 using 
infrared T B differences. T B difference techniques using MODIS will be very sensitive to 
thin clouds if the surface emittance and temperature are adequately characterized. Given 
a surface at 300°K and a cloud at 220°K, a cloud emittance of 0.01 increases the infrared 
window Tg by 0.5°K. Since the anticipated noise equivalent temperature for many 
MODIS infrared window spectral bands is better than 0.05°K, the cloud mask should be 
very sensitive to thin cloud (Ackerman et al., 1997). 

(1) Tri-Spectral Technique. Ackerman et al. (1997), suggests the 
use of a tri-spectral technique with T B (8pm) - T B (1 1pm) and T B (1 1pm) - T B (12pm) tests 
for thin cirrus clouds when precipitable water (pw) is greater than 1.0cm. A polar 
atmosphere is relatively dry and pw is less than 1.0cm. Therefore, only the T B (llpm) - 
T B (12pm) test is used to detect thin cirrus clouds. As will be discussed in the next 
chapter, MAS data suggests that this test might not be viable in daytime polar ocean 
regions. 
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(2) TB(3.9pm) - TB(ll|im). This test is used to identify fog and 
low clouds. Clouds at 11pm wavelength, independent of water droplet size, are very 
close to being black bodies, with emittance approximatley equal to 1. Clouds at 3.9pm 
wavelength are not black bodies and emittance is a function of water droplet size. For a 
normal water droplet size of 10pm, the emittance is approximately 0.85. 

During the day, the TsOlpm) is composed of only energy that has 
been emmitted from the cloud and is very close to the actual cloud top temperature. The 
TB(3.9pm) during the day is composed of both energy that is emmitted and reflected from 
the cloud. This will raise the sensed temperature above the actual cloud top temperature. 
Therefore, a positive difference in the TB(3.9pm) - Tb(1 1pm) test is an indication of 
liquid water clouds during the day. 

At night, negative values for TB(3.9pm) - TB(llpm) indicate a 
liquid water cloud. TB(3.9pm) is only a function of emmitted energy and will be less 
than the actual temperature since the emittance is less than one. Therefore, small 
T B (3.9pm) - T B (llpm) differences are observed when clouds are optically thick due to 
reduced emmitance at TB(3.9pm). 

c. Group 3 

The objective of Group 3 is to detect low clouds by using solar reflectance 
and reflectance ratio tests and is especially sensitive to thick, low clouds. It complements 
Group 1 because it is sensitive to thick, low level clouds while Group 1 has difficulty 
with low clouds that have small thermal contrast between cloud and background. This 
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technique works well when there is a high contrast in the reflectance between the surface 
and the cloud, for example, clouds over water. 

(1) Po.88|xm- This is a single band threshold test that does well for 
discriminating bright clouds over dark surfaces (e.g., stratus over ocean) and poorly for 
clouds over bright surfaces (e.g., sea ice). The reflectance test is adjusted for viewing 
angle and the test is eliminated when a sun glint condition exists for the pixel (Ackerman 
etal., 1997). 

(2) Ratio of po. 87 |un/po. 66 nm- Many earth surface features are less 
reflective at visible wavelengths than they are at near-IR wavelengths. Clouds do not 
exhibit these same variations in reflectance to any great degree (Ackerman et al., 1997). 

d. Group 4 

The objective of Group 4 is to detect high thin cirrus clouds during the 
day. Due to the fact that most of the atmosphere’s water vapor is located in the low to 
mid troposphere, most of the reflectance from the surface, low, and mid level clouds, 
should be attenuated prior to reaching any space borne sensor when sufficient water 
vapor (0.4 cm precipi table water) is in the atmosphere. Gao et al. (1993) developed a 
technique to detect the presence of thin cirrus clouds by a strong water vapor absorption 
band in the near-IR at 1.38|im (MAS - 1.88pm) wavelength. This band receives large 
amounts of scattered solar radiance when thin cirrus clouds are present. Ackerman et al. 
(1997) placed this technique in to its own group because they felt “there was no test that 
has been found to be as sensitive to the presence of thin cirrus.” 
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e. Group 5 

The objective of Group 5 is similar to those of Group 2 in that it is 
designed to detect thin clouds. However, in this case the focus is on thin cirrus clouds at 
night. In the polar ocean regions, a T B (3.9|i.m) - T B (1 1pm) difference is employed to 
detect cirrus clouds and is relatively insensitive to the amount of water vapor in the 
atmosphere (Hutchison and Hardy, 1995). This test is greater than zero in 
semitransparent cirrus as sub-pixel warm features dominate the shortwave window 
radiances within a pixel. 
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HI. MAS/MODIS CLOUD MASK 



MAS offers the opportunity to test the validity of the cloud mask prior to the 
MOD1S launch. Subjective analysis of the cloud properties are compared with cloud 
mask results. Modifications of tests and thresholds are then suggested to improve the 
cloud mask performance. 

A. OPTIMIZING THE CLOUD MASK 

MAS Flight #97-73 is a daytime flight from the FIRST ISSCP Regional 
Experiment - Arctic Cloud Experiment (FIRE/ACE) flight of June 2-3, 1998 and is used 
as a case for cloud mask optimization in the polar ocean region. Figure 3.1 shows the 
over-ground track of the ER-2 aircraft. 




Figure 3.1. The ER-2 over-ground track is plotted for MAS Flight #98-73, June 2-3, 
1998. (Gumley et al., 2000) The two regions of the flight used for cloud 
mask optimization are identified as blue squares. 
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1 . 



Original Cloud Mask 



Figure 3.2 is an annotated RGB image (Red: Band 20, Green: Band 10, Blue: 
Band 2) of track 17, center time of 01:04:29 UTC with Pt Barrow to the top of the image 
(South). The scene contains an area of mid to low clouds over a pack of mixed multi- 
year and first-year sea ice that dominates most of the scene. The Marginal Ice Zone 
(MIZ) is the transition area between pack ice and lesser concentrations of ice. Take note 
of the belting of sea ice jetting into water from the main pack ice. An area of 40 - 60% 
concentration of multi-year and first-year sea ice is located in the upper left hand comer 
of the scene. 

The original cloud mask classifies the scene as one of two types, Ocean Polar 
Snow Day or Ocean Day. This classification is very dependent on the Hall Test 
described in Chapter II. Table 3.1 lists the original group, test and thresholds for the two 
scene types. 



Table 3.1. Original Groups, Tests, and Thresholds 



Polar Ocean Snow Day 


Ocean Day 


Group 


Test 


Threshold 

(a,p,y) 


Group 


Test 


Threshold 

(a,p,Y) 


1 


T B (13.9pm) 


219,220,221 


1 


T B (ll|im) 


267, 270, 273 


T B (13.9pm) 


219,220,221 


2 


T B (3.9pm) - T b (1 1pm) 


11.0, 9.0, 7.0 


2 


T B (llpm)-T B (12pm) 


Dynamic 


T B (1 1pm) - T B (3.9gm) 


-10, -8, -6 








3 


Po.88^m 


8, 7, 6.5 


Po.87nm/p0.66^im 


0.95, 1.10, 1.15 
0.95, 0.90, 0.85 


4 


Pl.38iim 

(MAS - pi.88iim) 


4.0, 3.5, 3.0 


4 


Pl.38nm 

(MAS - Pl.88nm) 


4.0, 3.5, 3.0 
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Figure 3.2. 




MAS Flight #97-73, Track 17, Center Time of 01:04:29 UTC is displayed 
as an RGB image where Red is band 20, Green is band 10, Blue is band 2. 
Clouds, multi-year and first-year sea ice and partial ice concentration area 
are noted. 



Figure 3.3 shows the original cloud mask correctly classifying most of the clouds 
as Cloudy. It also correctly classifies the clear region over the pack ice as Confident 
Clear. The cloud mask encounters difficulty along the M1Z, where it classifies the clear 
areas as either Undecided or Cloudy. It also has trouble in the water areas where the 
pixel is clear of clouds and sea ice, but the mask still classifies it as either Probably 
Clear, Undecided or Cloudy. 

a. Hall Test 

The Hall Test is very important in determining which scene type is used in the 
cloud mask. Since the entire scene is below 280°K for TB(ll)im), the Hall Test is 
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initiated everywhere in the scene. Figure 3.4 shows the NDSI and po.xxpm tests for this 
scene. With the exception of a few cloudy pixels, the NDSI test classifies most of the 
scene as possibly being sea ice. The po.xxnm test checks the scene reflectance and 
classifies the polar water area as free of sea ice if the reflectance is less than 0.10 (green 
areas in Figure 3.4). 



Figure 3.3. 




□Confident Clear 
□Probably Clear 
^Undecided 
□Cloudy 



Original Cloud Mask is applied to MAS Flight #97-73, Track 17, Center 
Time of 01:04:29 UTC. Mask classifications are described in the legend. 



Figure 3.5 shows that the Hall Test correctly identifies the pack ice and 
much of the cloudy region. The test identifies some of the clouds as sea ice but as will be 
seen latter, are eventually classified as Cloudy when the final confidence value is 
computed. Finally, the floating sea ice in the partial ice concentration area is correctly 
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NDSI 



Figure 3.4. 




The NDSI test (a) and the po.88nm test (b) is applied to MAS Flight #97-73, 
Track 17, Center Time of 01:04:29 UTC. The figures are colored 
according to the thresholds set for each test. The NDSI test classifies 
pixels with values greater than 0.4 as Probably Ice. The po.88nm classifies 
pixels with values greater than 0.1 as Probably Ice. 
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identified as sea ice. Therefore, the Hall Test properly classifies pixels and initiates the 
groups and tests that need to be conducted for cloud detection. 




□Sea Ice 



Figure 3.5. MAS Flight #97-73, Track 17, Center Time of 01:04:29 UTC, Hall Test. 
Mask classifications are described in the legend. 



b. Individual Pixels 

Figure 3.3 reveals the greatest misidentification in this scene that occurred 
along the MIZ and the partial ice concentration area. Two representative pixels are 
selected to investigate this error. They are 1 and 2 pixels away, respectfully, from the 
main pack ice. Table 3.2 lists test results for a pixel that is classified as sea ice by the 
Hall Test and uses the groups and tests that correspond with Polar Ocean Snow Day. 
Table 3.3 lists test results for a pixel that is immediately adjacent to the first pixel and is 
classified as water and therefore uses the groups and tests that correspond with Ocean 
Day. 
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Table 3.2. 



Polar Ocean Snow Day Pixel Values for MAS Flight #97-73, Track 17, 



Scanline 


124304 


Sample 


305 


Group 


Test 


Threshold 

(a, p, y) 


Value 


Confidence 


1 


T B (13.9|im) 


219,220,221 


254.87 


1 


2 


T B (3.9|im) - T b ( 1 1 jxm) 


11.0, 9.0, 7.0 


-2.72 


1 


4 


Pl.88um 


4.0, 3.5, 3.0 


0.52 


1 


Overall Confidence (Q) 


1 



Table 3.3. Ocean Day Pixel Values for MAS Flight #97-73, Track 17, Center Time 
of 01:04:29 UTC, Scanline 124304, Sample 306 



Scanline 


124304 


Sample 


306 


Group 


Test 


Threshold 

(a, P, Y) 


Value 


Test 

Confidence 


Group 

Confidence 


1 


T B (1 lpm) 


267, 270, 273 


270.76 


0.63 


0.63 


T B (13.9nm) 


219,220,221 


255.13 


1 


2 


T B (llpm) - 
T B (12|im) 


1.136656, 

0.636656, 

-0.61334 


-0.05 


0.77 


0.77 


T B (3.9|im) - 
T B (llUm) 


-10, -8,-6 


2.29 


1 


3 


Po.88pm 


8, 7, 6.5 


7.93 


0.03 


0.03 


Po.87p.rn/ Po.66pm 


0.95, 1.10, 1.15 
0.95, 0.90, 0.85 


0.83 


1 


4 


Pl.88um 


4.0, 3.5, 3.0 


0.40 


1 


1 


Overall Confidence (Q) 


0.36 



When the cloud mask classifies the pixel as sea ice (Table 3.2), an overall 
confidence value of 1 is generated for the pixel. Even though the wrong scene type was 
initially selected, the cloud mask properly identifies the pixel as Confident Clear. When 
the cloud mask classifies the pixel as water (Table 3.3), the overall confidence value of 
0.36 is generated for the pixel. The correct scene type is selected but the cloud mask 
incorrectly identifies the pixel as Cloudy. The confidence values range from Possibly 
Clear to Cloudy when the surface is water and the confidence value is Confident Clear 
when the surface is sea ice. 
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The Hall Test may identify the pixel as water when a pixel possesses sub- 
pixel sea ice. The po.ssum test threshold for the Ocean Day scene type is sensitive enough 
to detect the increased reflectance from the pixel due to subpixel sea ice. If the 
reflectance is too great for water, the clear confidence level from the po.88nm test 
decreases. This occurs in the MIZ where sub-pixel sea ice occurrence is high. Tb( 1 1pm) 
and TB(llpm) - TB(12pm) tests also have low clear confidence levels (0.63 and 0.77 
respectfully) but are not as dependent on the presence of sub-pixel sea ice and are 
therefore better candidates for the threshold adjustments necessary for optimized cloud 
detection in polar ocean regions. 

c. T B (11/Jm) 

The TsOlpm) test produced a 0.63 clear confidence level that the pixel 
was clear for the Ocean Day scene type case. Figure 3.6(a) shows the RGB image of the 
scene with the Tb(1 1pm) test result overlaid where the test was conducted. The clouds 
are not cold enough to be classified as Cloudy according to the original threshold, while 
the partial ice concentration area is not warm enough to be classified as Confident Clear. 

Surface temperature controls the value of T B (1 1pm) when corrected for 
surface emittance and given low water vapor content. Hall et al. (1998) used data from 
MAS images collected in April 1995, near St. Lawrence Island, Alaska, and estimated the 
mean seawater temperature to be 271.6°K using seawater emittance of 0.985. Based on 
direct measurement, the freezing point of polar ocean water was found to be 271.4°K. 
Therefore, any pixel can be classified as possibly sea ice if its surface temperature is less 
than or equal to 271 ,4°K, the freezing point of polar seawater. 
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Figure 3.6(b) is similar to 3.6(a) except the NPS modified threshold for 
Tb( 1 l|xm) test is applied. The threshold was lowered by 1°K from the original threshold. 
Polar ocean water is now warm enough to be identified as clear (white), while clouds are 
still cold enough to be identified as clouds since they are colder that the (3-threshold ((3 = 
272°K). 

When the NPS threshold is applied to the TsOlM-m) test and the cloud 
mask is run, the clear confidence level of the test pixel was increased to 0.79 from the 
previous value of 0.63. Figure 3.7 shows the original cloud mask (a) compared with the 
cloud mask with the modified NPS Tb(1 1pm) threshold (b). The clear confidence level 
increases for many of the pixels in the partial ice concentration area. The clouds are still 
identified as Cloudy. However, an unacceptable number of pixels are still classified as 
Undecided. 

<L T B (ll/jm) - T B (12/jm) 

Figure 3.8 shows the Tb(1 1pm) - TB(12pm) test colorized for an average 
view angle and Tb(1 1pm) threshold. It is difficult to obtain any information about the 
scene from this test. Therefore, this test is removed from the Ocean Day scene type tests 
in the optimized cloud mask. 

Figure 3.9 shows the original cloud mask (a) compared with the cloud 
mask without the Tb(1 lpm) - T B (l2pm) test. There is significant improvement of the 
cloud mask in the partial ice concentration area showing fewer pixels as Cloudy or 
Undecided. Most of the water region has changed to Probably Clear with much of the 
noise is removed 
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— 1273 00 

A 



272.00 



271.00, 



Ongmal T B (11^m) Threshold 



NPS T Q (1lMm) Threshold 

Figure 3.6. Original Tb(1 1 f-im) test threshold (a) and the NPS Tb(1 lpm) test threshold 
(b) are applied to MAS Flight #97-73, Track 17, Center Time of 01:04:29 
UTC and overlaid on a RGB images. The figures are colored according to 
the thresholds set for each test. The original threshold classifies values 
greater than 270°K as Probably Clear. The NPS threshold classifies 
values greater than 269°K as Probably Clear. 
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□Confident Clear 
□Probably Clear 
BBUndecided 
nCloudy 



Cloud Mask T B (1lpm) Test. -Ik 

Figure 3.7. Original cloud mask (a) and the NPS cloud mask (b) are applied to MAS 
Flight #97-73, Track 17, Center Time of 01:04:29 MAS Flight #97-73, 
Track 17, Center Time of 01:04:29 UTC. Mask classifications are 
described in the legend. 



Original Cloud Mask 



□Confident Clear 
□Probably Clear 
^Undecided 
□Cloudy 
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Figure 3.8. 




The T B (ll|Xm) - TB(12fim) test is applied to MAS Flight #97-73, Track 
17, Center Time of 01:04:29 UTC. The scene is colorized for an average 
threshold. 
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Original Cloud Mask 



Figure 3.9. 




Belt of Ice from Mam Pack 



^□Confident Clear 
□Probably Clear 
□Undecided 
□Cloudy 




Cloud Mask without T 3 (11jtm) - T g (12^im) Test 



Original cloud mask (a) and the NPS cloud mask (b) are applied to MAS 
Flight #97-73, Track 17, Center Time of 01:04:29 MAS Flight #97-73, 
Track 17, Center Time of 01:04:29 UTC. Mask classifications are 
described in the legend. 
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T B (3.9/jm) - T B (ll/jm) 



e. 

Once the new threshold was defined for T B (llp.m) and the Tb(1 l|im) - 
T B (12|im) test was removed, other scenes indicated the T B (3.9p.m) - T B (lljim) test 
threshold might not be sensitive enough to detect low clouds. To be consistent 
throughout the algorithm, the T B (1 l(im) - T B (3.9|J.m) test in the Ocean Day scene type 
was changed to T B (3.9jim) - T B (1 ljim). The only difference in the test is a sign change. 

Figure 3.10 shows Track 9 of FIRE/ ACE Flight #97-73, with a center time 
of 23:00:22 UTC and is an example where the entire scene is classified as sea ice on the 
surface and is dominated by thick low to mid-level clouds above. Dark shadows and 
semi-transparent clouds near the top of the scene might affect the performance of the 
cloud mask. Finally, along the bottom, multi-year sea ice is visible through the thin, 
broken low-level clouds. 

Figure 3.11 shows the original cloud mask. The original cloud mask over- 
characterized the scene as Confident Clear, especially in the dark shadow and semi- 
transparent area where high clouds are casting a shadow over low clouds. The original 
cloud mask also had difficulty detecting the thin, broken, low-level clouds. 

The cloud mask modified with the new T B (1 lpm) threshold and with 
T B (ll{im) - T B (12|im) removed was run but there were no changes in the detection of 
clouds. This is due to the fact that the entire scene was classified as sea ice and only used 
Polar Snow Day as the scene type. The pi.gg and T B (13.9|i.m) thresholds correctly 
determined the scene contained no high clouds. The T B (3.9jim) - T B (1 1pm) test is 
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therefore the only candidate for threshold modification in order to correctly classify the 
low and mid-level clouds apparent in Figure 3.10. 




Figure 3.10. MAS FIRE/ACE Flight #97-73, Track 9, Center Time 23:00:22 UTC, is 
displayed as an RGB image where Red is band 20, Green is band 10, Blue 
is band 2. Thick low to mid-level clouds dominate the scene with dark 
shadows and semi-transparent clouds near the top. Finally, along the 
bottom, multi-year sea ice is visible through the thin, broken low-level 
clouds. 



The TB(3.9pm) - Tb(1 1pm) test has a flipped threshold. The test indicates 
the scene is Confident Clear until it reaches the lower threshold value, and the confidence 
of the pixel being clear decreases as the test reaches the upper threshold value. At this 
point, there is zero confidence that the pixel is clear. 
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Figure 3.12(a) shows the original T B (3.9|im) - T B (llpm) test threshold. 
The original threshold was unable to identify many of the low clouds in the scene. It also 
failed to identify the clouds in the dark shadow and semi-transparent area. 



Figure 3-11. 




□Confident Clear 
□Probably Clear 
SUndecided 
nCloudy 






Low - Mid Level Clouds, 
Semi-Transparent 



Original Cloud Mask is applied to MAS FIRE/ACE Flight #97-73, Track 
9, Center Time 23:00:22 UTC. Mask classifications are described in the 
legend. 



Figure 3.12(b) is similar to 3.12(a) except the NPS modified threshold for 
T B (3.9pm) - T b (1 1pm) test is applied. The new threshold reduced a from 1 1 to 5.25, (3 
from 7 to 3.5 and y from 7 to 1.75. The new threshold improved the detection of the low 
level clouds in the thin broken and dark shadow, semi-transparent areas. The NPS 
T B (3.9|lm) - T b (1 1pm) threshold still had trouble where the dark shadows reduced the 
difference enough to be falsely considered clear. 

Figure 3.13 shows the cloud mask with the NPS modified T B (3.9pm) - 
T B (1 lpm) test threshold added. The cloud mask detects more low and middle level 



36 




Original T B (3 9um) - T B (11nm) Threshold 






i ^ Low - Mid Level C 

Semi-Transparent ^ 

>v 

K i 

< • 



11.00 

▲ 



1 

1^-iL 




NPS T 8 (3 9j.tm) - T B (11jam) Threshold 



3.00 



2.50 



2.00 



5.00 



4 50 * 
3 
O 

O 

4.00 > 



Low - Mid Level Clouds, 
Semi-Transparent 



K *** • 



Figure 3.12. Original TB(3.9|im) - Tb(1 l|im) test threshold (a) and the NPS TB(3.9|im) 
- Tb(1 lpm) test threshold (b) are applied to MAS FIRE/ ACE Flight #97- 
73, Track 9, Center Time 23:00:22 UTC. The figures are colored 
according to the thresholds set for each test. The original threshold 
classifies values greater than 9.0 as Probably Cloudy. The NPS threshold 
classifies values greater than 3.5 as Probably Cloudy. 
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Figure 3.13. 




Low - Mid Level Clouds, 
j* Semi-Transparent 



.□Confident Clear 
□Probably Clear 
iUndecided 
□Cloudy 



Original Cloud Mask 




□Confident Clear 
□Probably Clear 
ESUndecided 
□Cloudy 




NPS Cloud Mask 



Original cloud mask (a) and the NPS cloud mask (b) are applied to MAS 
Flight #97-73, Track 9, Center Time of 23:00:22 UTC. Mask 
classifications are described in the legend. 
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clouds. Only the darkest shadow, semi-transparent area pixels are still classified as 
Confident Clear. 

B. NPS MAS/MODIS CLOUD MASK 

The observed weaknesses in the original cloud mask were revealed by two new 
scene type subroutines {Polar Ice Day and Polar Ocean Day). These scene types were 
developed after the removal of the T B (1 1pm) - T B (12pm) test and by identifying the 
limitations in the thresholds for T B (1 1pm), and T B (3.9|im) - T B (1 1pm) tests. These new 
scene types only run when the pixel is located in a daytime polar ocean region or when 
sea ice has been detected during the day. Table 3.4 lists the groups, tests and thresholds 
for these new scene type subroutines. 



Table 3.4. New Scene Type Subroutines. 



Polar Ice Day 


Polar Ocean Day 


Group 


Test 


Threshold 
(a, p, y) 


Group 


Test 


Threshold 
(a, p, Y) 


i 


T B (13.9|im) 


219,220,221 


1 


T B (llpm) 


266, 269, 272 


T B (13.9pm) 


219,220,221 


2 


T B (3.9|im) — T B ( 1 lpni) 


5.25,3.5 1.75. 


2 


T b (3 .9pm) — T b ( 1 1 pm) 


10, 8, 6 








3 


Po.88pm 


8, 7, 6.5 


Po.87pm/Po.66pm 


0.95, 1.10, 1.15 
0.95, 0.90, 0.85 


4 


Pi. 38pm 

(MAS - P 1.88pm) 


4.0, 3.5, 3.0 


4 


P 1.38pm 

(MAS ~ Pi. 88pm) 


4.0, 3.5, 3.0 



Figure 3.14 shows the final NPS cloud mask for Track 17, center time of 01:04:29 
UTC with the new scene type subroutines. The clouds are better defined and cover a 
larger area. A few pixels over the cloud-free pack ice are now misclassified as 
Undecided due to the change in the T B (3.9pm) - T B (1 1pm) test threshold and appears to 
be related to noisy pixels. 
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The new scene type subroutines greatly improved the classification of the partial 
ice concentration area. Most pixels are classified as either Confident Clear or Probably 
Clear - a direct result of the Tg(l lpm) threshold change and the removal of the 
Tb(1 lpm) - TB(12jim) test. 

The NPS cloud mask still has problems with the po.ssum test when pixel coverage 
is Polar Ocean Day and sub-pixel sea ice is present. Under these conditions, the 
reflectance value rises above the threshold, but changing this threshold degrades the 
ability of the Polar Ocean Day scene type to correctly classify the pixel as Confident 
Clear when no sub-pixel sea ice is present. Misclassification is most evident along the 
MIZ where there is a large concentration of sub-pixel sea ice. The artifact can be used to 
identify the MIZ or other sea ice edge features. Unfortunately, the problems associated 
with sub-pixel sea ice contamination will increase as the resolution of the pixel decreases 
from 50m (for MAS) to 1km (for MODIS). On the other hand, much of the noise in the 
MAS IR bands will decrease because the comparable MODIS IR bands are expected to 
have better radiometric sensitivity and more accurate calibration. 



40 



Figure 3.14. 
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Multi and First Year ice 



Original Cloud Mask 



□Confident Clear 
□Probably Clear 
ESUndecided 
□Cloudy 



Beit of Ice from Main Pack 



Original cloud mask (a) and the NPS cloud mask (b) are applied to MAS 
Flight #97-73, Track 17, Center Time of 01:04:29 UTC. Mask 
classifications are described in the legend. 
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IV. MAS DATA 



With the new thresholds in place, other MAS scenes were reviewed to test the 
optimization of the cloud mask. MAS data from the Winter Cloud Experiment 
(WINCE), First ISSCP Regional Experiment - Arctic Cloud Experiment (FIRE/ACE) and 
Winter Experiment (WINTEX) was used to verify the two new scene types. Subjective 
cloud analysis was done to verify the effectiveness of the new NPS MAS/MODIS cloud 
mask (NPS cloud mask). 

A. MAS FIRE/ACE FLIGHT #98-68, TRACK 6, CENTER TIME 21:04:36 
UTC 

The MAS FIRE/ACE Flight #98-68 (May 26-27, 1998), track 6, center time 
2 1 :04:36 UTC is shown in Figure 4. 1 . 




Figure 4.1. MAS FIRE/ ACE Flight #98-68, May 26-27, 1998, Flight Track. (Gumley 
2000) 
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This case is an example of a relatively clear arctic condition with a high solar 
zenith angle (57.9°). Figure 4.2 shows the scene is dominated by sea ice with thin clouds 
along the right hand side. The Hall Test classified the entire scene as sea ice, therefore 
only the Polar Ice Day scene type was used to detect cloud. The T B (13pm) test detected 
no high clouds in the scene. 



Figure 4.2. 




MAS FIRE/ ACE Flight #98-68, Track 6, Center Time 21:04:36, is 
displayed as an RGB image where Red is band 20, Green is band 10, Blue 
is band 2. Sea ice dominated scene with a thin cloud along the right hand 
side. 



Figure 4.3 shows the original cloud mask. It classifies most of the scene as clear 
with a cloud along the right side of the scene. The striping in the scene is directly related 
to the noise in the MAS IR bands. 
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Figure 4.3. 




Original Cloud Mask is applied to MAS FI RE/ACE Flight #98-68, Track 
6, Center Time 21:04:36 UTC. Mask classifications are described in the 
legend. 



Figure 4.4(a) shows the original T B (3.9|im) - T B (ll|lm) test threshold. The 
original threshold fails to detect a large portion of the thin cloud. Figure 4.4(b) shows the 
result of the NPS T B (3.9pm) - T B (1 l|im) test. The NPS threshold detects much more of 
the thin cloud and possibly a very thin cloud along the bottom. The NPS threshold 
results in significantly more noise in the clear area, where it has lowered the final clear 
confidence value in many of the clear pixels. Again, the lowered values are an artifact of 
the MAS 1R calibration. 

Figure 4.5 shows the original and NPS cloud masks. The NPS mask still 
classifies most of the scene as Confident Clear. The NPS threshold results in 
significantly Undecided pixels where it has lowered the final clear confidence value in 
many of the clear pixels. Identification of these Undecided pixels are an important result 
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Figure 4.4. 




Original T B <3.9pm) * T e (11nm) Threshold 




NPS T g (3 9^im) - T B (11nm) Threshold 



Original TB(3.9pm) - Tb( 1 ljurn) test threshold (a) and the NPS TB(3.9|im) 
- Tb( 1 1pm) test threshold (b) are applied to MAS FIRE/ACE Flight #98- 
68, Track 6, Center Time 21:04:36 UTC. The figures are colored 
according to the thresholds set for each test. The original threshold fails to 
detect a large portion of the thin cloud. The NPS threshold detects much 
more of the thin cloud and possibly a very thin cloud along the bottom. 
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Figure 4.5. 




Ongmai Cloud Mask 




NPS Cloud Mask 



□Confident Clear 
□Probably Clear 



□Cloudy 



Original cloud mask (a) and the NPS cloud mask (b) are applied to MAS 
FIRE/ACE Flight #98-68, Track 6, Center Time 21:04:36 UTC. Mask 
classifications are described in the legend. The NPS cloud mask classifies 
most of the scene as Confident Clear, and improves the detection of the 
thin cloud. 
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of the variable threshold techniques used in the MAS/MODIS cloud mask. Undecided 
pixels in Figure 4.5 identify regions of very thin clouds that are also apparent in Figure 
4.2. 

B. MAS FIRE/ACE FLIGHT #98-68, TRACK 16, CENTER TIME 23:31:38 
UTC 

Figure 4.6 shows the FIRE/ ACE Flight 98-68, track 16, center time 23:31:38 
UTC. This is an example of a scene that contains both sea ice and high clouds. Clouds 
are casting a shadow over the sea ice that, at lower spatial resolution, could be 
misidentified as a lead. As the main cloud thins along its edges, it becomes increasing 
difficult to distinguish between cloud and sea ice. 

Figure 4.7 shows the original cloud mask. It identifies the main cloud mass but 
has difficulty along the edge. There is no indication of the thin clouds in the center of the 
scene suggested in Figure 4.6. The Hall Test classifies the entire scene as sea ice, 
therefore only the Polar Ice Day coverage subroutine is run. Again, TB(13.9|im) did not 
detect any clouds in the scene. 

Figure 4.8(a) shows the original TB(3.9|im) - Tb(1 l(im) test threshold. The 
original threshold detects much of the main cloud mass but subjective analysis suggests 
that this test is the reason why the original cloud mask is having difficulties with the thin 
cloud edges. Figure 4.8(b) shows the NPS TB(3.9|im) - TB(ll|im) test threshold. The 
NPS threshold is detecting the thin clouds that were missed by the original threshold. It 
has also detected two smaller clouds along the bottom right of the scene. Although this is 
a daytime case and these two clouds are in shadow, the emittance difference between 
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Figure 4.6. 



Figure 4.7. 




MAS FI RE/ ACE Flight #98-68, Track 16, Center Time 23:31:38 UTC, is 
displayed as an RGB image where Red is band 20, Green is band 10, Blue 
is band 2. The scene contains both sea ice and high clouds, with clouds 
casting a shadow over the sea ice. 




□Confident Clear 
□Probably Clear 
SlUndecided 
□Cloudy 



Original Cloud Mask is applied to MAS FIRE/ACE Flight #98-68, Track 
16, Center Time 23:31:38 UTC, Mask classifications are described in the 
legend. 
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Tb(3.9(hti) and TB(llpm) still reveals the presence of the clouds. Again, noise is 
introduced into the test by the MAS IR bands. 

Figure 4.9 shows the NPS cloud mask for this case. It detects the thin clouds 
along the edge of main cloud mass due directly to the threshold modifications in the 
TB(3.9|im) - Tb(1 lpm) test. It also detects the two smaller clouds along the bottom right. 
Noise in the cloud mask is generated by the TB(3.9pm) - Tb(1 lpm) test. 

C. MAS FIRE/ACE FLIGHT #98-68, TRACK 4, CENTER TIME 20:11:05 

UTC 

MAS FIRE/ ACE Flight #98-68, Track 4, Center Time 20:11:05 UTC, is a 
difficult case because sea ice is visible through semi-transparent clouds. The case is 
shown in Figure 4.10 where an optically thick cloud is casting a shadow on the clouds 
below it. Gravity waves are detectable in the thin, low-level clouds. Tb( 13.9pm) 
detected very little cloud and was not a factor in the ultimate confidence determination. 

Figure 4.11 shows the original cloud mask. It classifies much of the scene as 
Cloudy or Undecided. It also classifies a portion of the shadow and the lower right 
comer of the scene as Confident Clear. 
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Figure 4.8. 




Original T 0 (3.9j.im) - T B (11pm) Threshold 




NPS T g (3 9fim) - T B (11pm) Threshold 



Original T B (3.9lim) - T B (1 l|im) test threshold (a) and the NPS T B (3.9|lm) 
- T b (1 lpm) test threshold (b) are applied to MAS FIRE/ACE Flight #98- 
68, Track 16, Center Time 23:31:38 UTC. The figures are colored 
according to the thresholds set for each test. The original threshold did not 
detect the thin clouds, while the NPS threshold improved the thin cloud 
detection. 
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Figure 4.9. 




Original Cloud Mask 




Clear 



Cloud 



Shadow 



Cloud 



NPS Cloud Mask 



□Confident Clear 
□Probably Clear 
^Undecided 
□Cloudy 



Original cloud mask (a) and the NPS cloud mask (b) are applied to MAS 
FIRE/ACE Flight #98-68, Track 16, Center Time 23:31:38 UTC. Mask 
classifications are described in the legend, (b) NPS Cloud Mask. The 
NPS cloud mask improved the detection of thin clouds. 
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Figure 4. 10. 




MAS FI RE/ ACE Flight #98-68, Track 4, Center Time 20:11:05 UTC, is 
displayed as an RGB image where Red is band 20, Green is band 10, Blue 
is band 2. Optically thick cloud is casting a shadow on the clouds below 
it. Sea ice is visible through semi-transparent clouds and gravity waves are 
detectable in the thin, low-level clouds 
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□Confident Clear 
□Probably Clear 
[HUndecided 
□ Cloudy 



Figure 4. 11. Original Cloud Mask is applied to MAS FIRE/ACE Flight #98-68, Track 
4, Center Time 20:1 1:05 UTC. Mask classifications are described in the 
legend. Much of the scene as classified as Cloudy or Undecided. A 
portion of the shadow as classified Confident Clear. The lower right 
comer of the scene, the cloud mask identifies the area as Confident Clear. 




Figure 4.12 shows where the Hall Test identifies sea ice. The pack ice in the 
scene and much of the multiyear ice floating in the open water is classified as sea ice. 
The test did not classify much of the optically thick cloud as sea ice. 

Figure 4.13 shows the TB(3.9pm) - Tb(11!Uii) difference. Semi-transparent low- 
level clouds and gravity waves are resolved. The test also resolves the shadow from the 
thick liquid water cloud. 

When the Hall Test classifies the pixel as sea ice, the Polar Ice Day coverage is 
used. Figure 4.14(a) depicts the original TB(3.9|0.m) - Tb(1 l|im) test threshold. The test 
classifies the shadow area and much of the low level clouds as clear. The original test 
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resolves portions of the gravity waves in the low level clouds but the differences do not 
exceed the threshold necessary to classify the entire region as cloudy. 




Figure 4.12. MAS FIRE/ACE Flight #98-68, Track 4, Center Time 20:1 1 :05 UTC, Hall 
Test overlaid on a RGB image. The pack ice in the scene and much of the 
multiyear ice floating in the open water is classified as sea ice. The test 
did not classify much of the optically thick cloud as sea ice. 



Figure 4.14(b) shows the NPS T B (3.9p.m) - T B (ll|im) test threshold. The test 
now correctly classifies the shadow area as cloud. Much of the low level clouds that 
were missed by the original threshold are now identified as cloud. Finally, the gravity 
wave features are fully resolved by the test. 
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Figure 4.13. T B (3.9pm) - Tb(1 1pm) difference applied to MAS FIRE/ACE Flight #98- 

68, Track 4, Center Time 20:1 1:05 UTC. The color scale on the right is 
the range of temperatures found in this scene. Semi-transparent low-level 
clouds and gravity waves are resolved. The test also resolves the shadow 
from the thick liquid water cloud. 



Another Polar Ocean Day test is TB(llpm). Figure 4.15(a) shows the original 
T B (llpm) test threshold. The test correctly classifies the scene as cloudy. It has no 
trouble in the shadow area classifying it as cloud. Figure 4.15(b) shows the NPS 
TB(lljim) test threshold. It raises the clear confidence level where thin clouds are 
present over the sea ice but not to a level where it might be classified as Confident Clear. 

Figure 4.16 shows the original and NPS cloud mask. The NPS cloud mask 
correctly classifies the shadow due to the threshold modifications in the TB(llpm) and 
T B (3.9pm) - T b (1 1pm) tests. The low-level clouds are also correctly classified as Cloudy 
due to the modification in the TB(3.9pm) - Tb(1 1pm) test. 
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Cloud 



NIPS T B {3.9pm) - T 0 (11prn) Threshold 

Figure 4.14. Original T B (3.9pm) - T B (1 lpm) test threshold (a) and the NPS T B (3.9pm) 
- T B (1 1pm) test threshold (b) are applied to MAS FIRE/ACE Flight #98- 
68, Track 4, Center Time 20:11:05 UTC and overlaid on T B (3.9|im) - 
T b (1 1pm) difference. The figures are colored according to the thresholds 
set for each test. The original threshold classifies the shadow area and 
much of the low-level clouds as clear. The NPS threshold correctly 
classifies the shadow and low-level clouds as clouds 
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Figure 4.15. 
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Original TB(llpm) test threshold (a) and the NPS (lljim) test threshold 
(b) are applied to MAS FIRE/ACE Flight #98-68, Track 4, Center Time 
20:11:05 UTC and overlaid on TeCllpm). The figures are colored 
according to the thresholds set for each test.. Both tests classify the scene 
as generally cloudy. 
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Figure 4.16. Original cloud mask (a) and the NPS cloud mask (b) are applied to MAS 
FIRE/ ACE Flight #98-68, Track 4, Center Time 20:11:05 UTC. Mask 
classifications are described in the legend. The NPS cloud mask correctly 
classifies the shadow and low-level clouds as Cloudy. 
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D. MAS WINCE FLIGHT #97-45, TRACK 6, CENTER TIME 20:18:05 UTC 

Figure 4.17 shows the ground track for Flight #97-45 conducted on February 6, 
1997. This flight in not in a polar ocean region but has lake ice and is in the area of 
responsibility of the NIC. 




FLIGHT 97-048 S FEBRUART 1997 R/C 708 MASSO 

LAMBERT CONFORMAL PROJECTION! SP1 = 40-7 SP2 = 53-0 CM = -84.4 ROTATED BT 0-0 

17*21 *14 TO 22 « 42 » 38 UT SCALE !«9.S3E-06 TIME TICK EVERT 10-00 MINUTES 



Figure 4. 17. MAS WINCE Flight #97-45, Flight Track. (Gumley, 2000) 

Figure 4.18 shows an RGB image of MAS WINCE Flight #97-45, track 6, with a 
center time of 20:18:05 UTC. It is in James Bay between the coast and Akimiski Island. 
There is fast ice along the shoreline with a large floe in the center of the scene. Leads are 
located between the fast ice and the floe where no sea ice is apparent. Small clouds are 
casting elongated shadows on the sea ice due to the high solar zenith angle of 75.2°. 

The Hall Test classified the entire scene as either sea ice or snow with the 
exception of small areas that correspond to visible clouds. This included the leads that 
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appear to be open water in the visible and the 1R bands. Figure 4.19(a) shows the NDSI 
test. With the exception of a few small areas where clouds are present, the NDSI value of 
the entire scene is greater than 0.4 and is therefore possibly sea ice or snow over land. 



Figure 4.18 




MAS WINCE Flight #97-45, Track 6, Center Time 20:18:05 UTC is 
displayed as an RGB image where Red is band 20, Green is band 10, Blue 
is band 2. Fast ice located along the shoreline with a large floe in the 
center of the scene. Leads are located between the fast ice and the floe 
where no sea ice is apparent. Small clouds are casting elongated shadows 
on the sea ice 



Figure 4.19(b) shows po.88nm f° r this scene. The entire scene may contain either 

snow or sea ice because no po.88fim value was less than 0.1. In the leads, where po.88|im is 

expected to be less than 0.1, there were no values less than 0.12. This may be due to a 

thin sheet of ice in the leads that is thin enough to be transparent at visible wavelengths 

but increases po.88pm enough to pass the test for sea ice. 
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Figure 4.20 shows the original cloud mask. It correctly identifies much of the 
scene as Confident Clear. The clouds were detected, but a subjective analysis suggests 
that the cloud mask could do a better job in depicting the size and shape of the clouds. 
With the lead classified as sea ice, the cloud mask correctly identifies it as Confident 
Clear. 

Figure 4.21 shows the TB(3.9pm) - T B (1 lpm) difference. This reveals the low 
level clouds in the scene. Figure 4.22(a) shows the original threshold as applied to the 
T B (3.9pm) - TB(lljJ-m) test. The original threshold detects some of the clouds. Figure 
4.22(b) shows the NPS threshold applied to the T B (3.9|im) - T B (1 1pm) test . Subjective 
analysis reveals that the NPS Ts(3.9pm) - Tb(1 1pm) test threshold increases the area of 
cloud detection. The clouds are better defined and better match the shadows in Figure 
4.18. When fast ice and sea ice are present, there is a decrease in the clear confidence 
level for this test. In leads, the test is negative, indicating that T B (1 1pm) is greater than 
T B (3.9pm), thus increasing the clear confidence level. Again, the MAS IR channels are 
noisy. 

Figure 4.23 shows the original and the NPS cloud mask. The new threshold for 
the T B (3.9pm) - TB(llpm) test has increased the number of pixels classified as Cloudy. 
Many pixels are classified as Undecided due to the noise in the TB(3.9pm) - TsOlpm) 
test, but these pixels still generate a confidence level of greater then 0.85 clear confidence 
level. 
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NDSI 



Figure 4.19 




The NDSI test (a) and the po.ssum test (b) is applied to MAS WINCE Flight 
#97-45, Track 6, Center Time 20:18:05 UTC. The figures are colored 
according to the thresholds set for each test. The NDSI test classifies 
pixels with values greater than 0.4 as Probably Ice. The Po ssum classifies 
pixels with values greater than 0. 1 as Probably Ice. 
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Figure 4.20 Original Cloud Mask is applied to MAS WINCE Flight #97-45, Track 6, 
Center Time 20:18:05 UTC. Mask classifications are described in the 
legend. 




Figure 4.21 T B (3.9|im) - T B (1 1pm) difference applied to MAS WINCE Flight #97-45, 

Track 6, Center Time 20:18:05 UTC. The color scale on the right is the 
range of temperatures found in this scene. 
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E. 



STATISTICAL ANALYSIS 



Table 4.1 shows the results of analyzing 18 tracks of data from three MAS flights. 



The average length of these tracks was 350 km and covered and average area of about 
13,000 km 2 . The NPS cloud mask classified 12.1% more pixels as Cloudy than the 
original mask. Subjective analysis of these tracks confirms that the new Cloudy pixels 
were indeed consistent with cloud features and not due to erroneous classification. 
Confident Clear pixels decreased by 35.2% in the NPS cloud mask indicating far fewer 



pixels are incorrectly identified as cloud-free. The NPS cloud mask generated about 
twice as many Probably Clear and Undecided pixels as the original cloud mask due to 



greater sensitivity to small, thin cloud features. 



Table 4.1 



Statistica 



summary of original and NPS cloud 



mask pixel classific ation 





Original 
Cloud Mask 


NPS 

Cloud Mask 


Change 


Total Pixels 


92388744 


92388744 




% Confident Clear 


39781316(43.0%) 


25769650 (27.9%) 


-14011666 (-35.2%) 


% Probably Clear 


1276156(1.4%) 


2444210 (2.6%) 


1168054(91.5%) 


% Undecided 


4039600 (4.4%) 


9200997 (10.0%) 


5161397 (127.8%) 


% Cloudy 


47291672 (51.2%) 


54973887 (59.5%) 


5512242(12.1%) 
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Original T e (3.9um) - T 0 (1^m) Threshold 



Figure 4.22 




NPS T B (3.9pm) - T e (11pm) Threshold 



Original T B (3.9pm) - T B (1 ljim) test threshold (a) and the NPS T B (3.9|im) 
- T B (1 1pm) test threshold (b) are applied to MAS WINCE Flight #97-45, 
Track 6, Center Time 20:18:05 UTC. The original threshold detects some 
of the clouds. The NPS threshold increases the area of cloud detection and 
they are better defined 
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Figure 4.23. 




nConfident Clear 



Original Cloud Mask 




□Confident Clear 
□Probably Clear 
□Undecided 
□Cloudy 



NPS Cloud Mask 



Original cloud mask (a) and the NPS cloud mask (b) are applied to MAS 
WINCE Flight #97-45, Track 6, Center Time 20:18:05 UTC. Mask 
classifications are described in the legend. 
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V. MODIS DATA 



Figure 5.1 shows Terra's eastern Arctic ground track for April 1, 2000. Granule 
17:00 UTC is a descending, daytime, track and is used to test the two new scene types on 
MODIS data. 




Figure 5.1 . Terra’s eastern Arctic ground tacks for April 1, 2000 (SSEC, 2000). 

A. CALIBRATION ISSUES 

There are calibration problems with some of the bands on the MODIS instrument 
and are believed to be correctable. Calibration problems with only the bands that affect 
the two new scene types are discussed below. 

1. Band 26 (1. 38pm) 

The objective of Band 26 (1.38pm) is to detect high cirrus clouds. Unfortunately, 
there is electronic cross talk with Band 26 and bands that see the surface. This allows 
surface emittance to contaminate the reflectance value for Band 26. 
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2 . 



Band 31 (1 1pm) 



Band 31 (11pm) lias optical cross talk contamination and is on the same focal 
plane as Bands 32 - 36. After 1 lp.ni energy passes through the filter, some of the energy 
does not reach the detector but remains in the substrate and contaminates the remaining 
Bands. 

3. Band 35 (13.9pm) 

Band 35 (13.9pm) has a scan mirror angle problem. Figure 5.2 is a graphic 
illustrating the scan mirror contamination. As the scan moves from left to right, there is 
an unexpected increase in the scan mirror angle. This leads to emittance from the scan 
mirror that increases the radiance value for Band 35. Therefore, an error is induced in 
T b (1 3.9pm) as the scan mirror tracks from left to right from a correct value to a warmer 
value. 

Scan Track 



Too Warm 
^ Granule 



Figure 5.2. Band 35 (13.9pm) Scan Mirror Error. As the scan mirror tracks left to 
right, the mirror angle increases, which increases the brightness 
temperature across the scan. 




Along Track 
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B. 



HUDSON BAY 



Figure 5.3 shows the April 1, 2000, 17:00 UTC, MODIS granule of the Hudson 
Bay. Figure 5.3(a) shows Band 2 (0.865fim) and Figure 5.3(b) shows Band 31(1 l|im). 
Hudson Bay is mainly cloud free with first year sea ice dominating the entire bay. Mid to 
low-level clouds are present along the west coast of the bay and extending northeastward 
between Coats and Mansel Islands. Finally, high clouds cover the Foxe Basin to the 
north and James Bay to the south. 

Figure 5.4 shows the TB(lllim) test for the April 1, 2000, 17:00, UTC MODIS 
granule with a land mask applied to highlight the water areas. Figure 5.4(a) shows the 
original threshold, while Figure 5.4(b) shows the NPS threshold. As described in Chapter 
III, the original MODIS Tb(1 lp.m) test threshold is too high for a polar ocean region. 
When the NPS threshold is applied to the MODIS TB(ll|im) test, the clear confidence 
level for an individual pixel became more representative for a polar ocean scene. 

At the location indicated by A in Figure 5.4, the NPS threshold improves the 
detection of the clear skies in the Davis Strait because the new threshold takes into 
account the lower surface temperatures associated with polar oceans. At B in Figure 5.4, 
the NPS threshold raises the clear confidence level over the sea ice in the Hudson Strait 
where some low and high clouds are present. Finally, at location C in Figure 5.4 the NPS 
threshold raises the clear confidence level over the sea ice where the low clouds are 
present. The increase in the clear confidence level at locations B and C is not enough to 
incorrectly identify the locations as clear. 
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Figure 5.5 shows the TB(3.9|nm) - TB(ll|JJn) test for April 1, 2000, 17:00 UTC, 
MODIS granule. Again, Figure 5.5(a) shows the original threshold and Figure 5.5(b) 
shows the NPS threshold. As with the MODIS T B (lllim) test, the MODIS Tb(3.9|Hii) - 
Tb(1 lpm) test threshold is too high for a polar ocean region. Locations A and D are 
regions where high clouds cover the surface. The NPS threshold lowered the clear 
confidence level for these regions enough for the high clouds to be classified correctly. 
For low clouds, locations B and C had the greatest improvements. The new threshold 
detected the low level cloud that was missed at B and increased the detection of low 
clouds at C. 

Figure 5.6(a) shows the original cloud mask. It classifies portions of the clear 
area in the Davis Strait (A) as Undecided to Confident Clear. Figure 5.6(b) shows the 
NPS cloud mask correctly classifies location A as Confident Clear. This is a direct result 
of lowering the T B (ll|Lim) test threshold. The NPS cloud mask also changed the 
classification of the clouds at location B so that almost all of the pixels are classified as 
Cloudy. At location C, where high clouds covered James Bay, the NPS cloud mask 
correctly increased the number of pixels classified as Cloudy. The changes at B and C 
were a direct result of lowering the TB(3.9jim) - T B (1 l|J.m) test threshold. The 
demarcation between Undecided and Confident Clear in the center of both cloud masks 
(under the label “Foxe Basin” and “Hudson Bay”) is due to Band 35 scan mirror error. 
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free with first year sea ice dominating the entire bay. Mid to low-clouds are present long the west coast of the bay and 
extending north eastward between Coats and Mansel Islands. Finally, high clouds cover the Foxe Basin to the north 
and James Bay to the south 
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with P as midpoint. Brown indicates where the pixel’s ecosystem was either coast or land and the NPS cloud mask as 
not applied 
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granule. Mask classifications are described in the legend. 




CONCLUSIONS AND RECOMMENDATIONS 



A. CONCLUSIONS 

Initial optimization of the MODIS cloud mask for the polar ocean regions is 
described based on 8 daytime MAS flights. 18 tracks (approximately 350km x 36km) 
were selected from 3 of these flights that had the greatest scene variability. MAS flights 
revealed that daytime thresholds for T B (lljim) and TB(3.9|im) - Tb(1 lpm) tests are too 
high for the polar ocean regions and must be lowered to compensated for lower average 
surface temperatures. In addition, the TB(ll|J.m) - T B (12|im) test gives little useful 
information during the daytime. Finally, Tb( 13.9pm) test was never a factor in the cases 
examined in this study. 

Because of the unique nature of the polar ocean regions, two new scene types 
were developed for the MODIS cloud mask. The new scene types incorporate new 
thresholds appropriate for polar ocean regions and are included in the NPS modified 
MAS/MODIS cloud mask (NPS cloud mask). The first new scene type is Polar Ice Day 
and is initiated during the day and when the pixel is identified as having sea ice. The 
second scene type is Polar Ocean Day and is initiated during the day and when the pixel 
is identified as water located greater than 60° north and south. 

The modifications greatly improve the detection of clouds over cold polar oceans 
where sub-pixel sea ice may be present or water temperatures might falsely indicate 
clouds. The number of Cloudy pixels (< 0.66 clear confidence level) for a given scene 
was increased 12.1% on average for MAS cases. The NPS cloud mask also classified 
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two times more Probably Clear and Undecided pixels than the original mask due to 
greater sensitivity to thin, small clouds. 

The improvements that were observed in the MAS data were confirmed when the 
NPS cloud mask was applied to 5 uncorrected MODIS cases. One case is presented here. 
The NPS cloud mask improved the clear confidence level when surface temperatures 
were below 273°K. It also lowered the clear confidence level when there were high and 
low clouds over the sea ice. 

The NPS MODIS cloud mask that was derived from Ackerman et al. (1997) will 
provide the NIC with the means to assist its analysts by giving them a tool to correctly 
identify clouds that might have been misidentified as sea ice in the past. A better quality 
sea ice analysis will result. Ultimately, this cloud mask will also lay the foundation for 
the development of ice recognition algorithms. 

B. RECOMMENDATIONS 

The operational MODIS cloud mask must have new scene types that are 
optimized for the polar ocean regions. It is not appropriate for thresholds, set for the 
globe, to be applied to the polar oceans due to their distinctly different surface 
temperatures. Unfortunately, further work should not proceed until calibration is 
completed for the MODIS instrument. Once completed, the two new scene types must be 
tested on a larger daytime dataset consisting of Arctic and Antarctic granules. Finally, 
work needs to be done on thresholds for nighttime cases that were not possible for this 
thesis. 
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